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Abstract 


The  U.S. “French  cooperative  work  in  the  field  of  Environmental 
Engineering  initiated  during  and  after  the  visit  of  the  author  to  the 
Centre  Scientifique  et  Technique  du  Batiment  (CSTB)  is  reported.  The 
topics  covered  in  this  report  are : 

1.  Comparative  results  of  digital  and  analog  calculations 
for  heating /coo ling  loads 

2.  Derivation  of  room  weighting  factors 

3.  In-depth  observations  of  activities  in  the  Hygrothermic 
Section  of  CSTB 

4.  Supplementary  visits  to  other  European  laboratories 


for  Environmental  Engineering. 


1.  Introduction 


During  the  visit  of  the  first  French  team  for  the  U.S. /French 
Cooperative  Program  headed  by  Mr.  Blachere  to  the  Building  Research 
Division  of  the  National  Bureau  of  Standards  in  November  of  1969,  it 
became  apparent  that  the  Environmental  Engineering  Section  of  NBS  and 
the  Hygrothermic  Section  of  CSTB  have  each  been  engaged  in  several  re- 
search programs  such  as  heating /coo ling  load  calculation,  ventilation, 
wall  heat  conduction,  and  room  air  distributions.  Both  organizations 
have  active  programs  for  improving  the  calculations  of  the  heating/ 
cooling  load  of  buildings.  While  the  effort  at  NBS  has  been  centered 
around  the  use  of  digital  computers,  French  counterparts  at  CSTB  liave 
been  employing  an  elaborate  analog  computer. 

Mr.  Borel,  who  represented  the  Hygrothermic  Section  of  CSTB,  and 
Mr.  F.  J.  Powell  of  the  Environmental  Engineering  Section  of  NBS, 
mutually  agreed  that  comparative  studies  of  heating/cooling  load  cal- 
culations for  the  same  building  performed  by  these  two  schemes  should 
be  of  great  interest  for  both  organizations. 
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In  addition,  it  was  felt  that  a researcher  from  the  Environmental 


Engineering  Section  would  benefit  very  much  by  observing  many  of  the 
French  research  activities  related  to  the  thermal  performance  of  build- 
ings. Subsequently,  the  author  was  selected  to  visit  the  CSTB  to  im- 
plement this  phase  of  the  U.S. /French  Cooperative  Program.  The  program 
was  carried  out  during  the  author’s  stay  at  CSTB  in  the  period  covering 
July  3 through  July  28  according  to  the  program  shown  in  Appendix  A. 

The  comparison  of  the  heating /coo ling  load  calculation  was  completed 
and  it  was  found  that  the  results  obtained  by  two  different  computing 
schemes  agreed  very  well.  Based  upon  this  experience  of  comparing  the 
computations  by  these  two  different  types  of  computers,  CSTB  engineers 
prepared  a new  definition  for  extending  the  cooperative  program  of  this 
kind,  which  is  shown  in  Appendix  B.  Also  included  in  this  report  are 
the  summaries  for  the  visits  to  several  other  European  laboratories. 
These  visits  took  place  during  and  after  the  author’s  stay  at  CSTB. 
Other  computed  results  of  the  previous  studies  made  at  the  CSTB  were 
utilized  to  derive  weighting  factors  needed  for  the  improvement  of  the 
NBS  digital  computer  program,  detail  of  which  is  shown  in  Appendix  C. 

At  the  end  of  the  visit  a proposal  stressing  the  need  for  a continued 

exchange  of  personnel  and  technical  information  was  drafted.  The 
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proposal  is  attached  to  Appendix  D. 
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2.  Migliiights  of  Llie  CSTB  Visit 


On  July  3,  1970,  the  author  arrived  in  Paris  and  started  to  work 
at  the  research  laboratories  of  CSTB  (Centre  Scientifique  et  Teclinique 
du  Batiment),  which  is  located  in  Champs -sur-Marne , approximately  20 
miles  east  of  Paris.  He  stayed  there  until  July  28  and  engaged  in 
technical  programs  arranged  by  the  CSTB  personnel.  The  chronological 
detail  of  the  program  is  attached  in  Appendix  A. 

During  the  author’s  stay,  he  worked  with  most  of  the  engineers  in 
the  Hygrothemic  Division  (Mr,  Croiset  - Head),  which  is  divided  into 
two  sections;  Thermal  Calculation  Section  (Mr.  Anquez  - Chief)  and  De- 
sign Concept  Section  (Mr.  Borel  - Chief). 

Following  are  the  highlights  of  his  program: 

1.  The  major  purpose  of  the  visit  was  to  check  the  thermal 
load  calculation  procedures  employed  by  CSTB  engineers 
on  their  analog  computer.  A fictitious  building  of  11 
rooms  located  at  the  latitude  of  42  %/  and  longitude  of 
88  °N  was  chosen  to  compare  the  cooling  load  under  a 
selected  weather  condition.  The  detail  of  building  con- 
struction and  occupancy  and  energy  usage  schedules  are 
all  shown  in  Figures  1 and  2 and  Tables  1>  2 and  3.  Only 
the  cooling  loads  of  5 rooms  on  the  second  floor  of  the 
building  were  considered  for  the  comparison.  Figures 
3 and  4 show  the  daily  profiles  of  the  cooling  load  of 
room  11  obtained  by  the  U,  S.  Post  Office  digital  computer 
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program  and  that  calculated  by  the  CSTB  analog  computer 
respectively,  which  shows  very  good  agreement  for  the 
maximum  thermal  loado 

Apparent  discrepancy  in  the  noon  time  dips  of  the  thermal 
load  of  the  two  calculations  is  due  to  the  fact  that  the 
internal  heat  load  was  treated  as  instantaneous  in  the 
CSTB  calculations  whereas  the  weighting  factors  to  distribute 
it  in  a time  domain  were  employed  in  the  NBS  calculations. 
Similar  good  agreements  between  the  digital  and  analog  com- 
putations were  obtained  for  the  other  four  rooms  of  the  same 
sample  building.  In  addition  to  the  comparison  of  the 
thermal  load,  room  temperature  calculations  for  the  NBS  ex- 
perimental building  were  also  made  by  the  NBS  digital  com- 
puter and  CSTB  analog  computer.  The  agreement  for  the  room 
temperature  calculation  was  perfect.  The  CSTB  analog  com- 
putation is  based  upon  a very  exact  radiative  heat  exchange 
process  among  the  interior  surfaces  of  the  room  whereas  the 
U.  S.  Post  Office  computer  program  approximates  this  process 
by  the  weighting  factor  calculation  procedure.  The  good 
agreement  obtained  between  the  two  calculations  is  ex- 
tremely important  for  NBS,  who  was  evaluating  the  U.  S. 

Post  Office  program  under  a specific  research  contract, 
with  the  U,  S,  Post  Office  Department. 
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2.  Mto  Moyer  has  excellent  background  in  analyzing  heat  con- 
duction problems  with  the  use  of  electrolytic  tank  and 
electric  conduction  paper.  These  devices  were  being  used 
to  study  heat  conduction  from  underground  tunnels,  exposed 
building  floors  (cold  bridge  effect),  and  underground  heat 
distribution  systems. 

3.  M.  Gibert*s  activities  on  the  ventilation  of  flats 
(apartment  houses)  and  air  flow  characteristics  of 
various  air  supply  and  exhaust  systems  are  parallel 
to  the  research  activities  initiated  in  NBS  by  Mr. 
Myklebost  and  being  continued  by  Professor  Hill. 

Although  the  basic  aims  of  these  two  investigations 
are  not  identical,  it  is  possible  to  exchange  ex- 
perimental informationo  The  goal  of  the  NBS  study 
is  to  obtain  dynamic  characteristics  of  the  air 
leakage  in  conjunction  with  the  outdoor  wind  charac- 
teristics and  building  construction  (crack  length  and 
width),  while  M.  Gibert's  objective  was  to  determine 
average  air  leakage  rate  from  flats  as  functions  of 
average  outdoor  weather  condition  and  window  openings. 

Mr.  Gibert  operates  an  experimental  flat  at  Veyre,  in 
the  vicinity  of  CSTB  research  laboratory.  An  Agfa 
motorized  camera  was  used  to  observe  the  window  opening 
habits  of  the  residents  in  the  apartment  house. 
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4.  An  optimization  study  on  the  heat  recovery  system  bein:^ 
conducted  by  M.  Galivier,  one  of  the  digital  computer 
oriented  engineers  at  CSTB,  was  very  interesting.  He 
proposes  an  unique  heat  exchanger  for  the  ventilation 
and  exhaust  air  streams  which  utilizes  multiple  layers 
of  corrugated  plastic  sheets  in  a cross-flow  fashion. 

The  size  of  the  heat  exchanger  was  determined  so  that 
heat  transfer  effectiveness  and  fan  horsepower  require- 
ment can  be  optimized. 

5.  Mr.  Bertolo  has  developed  a printed  circuit  long-wave- 
length radiometer  to  measure  the  radiative  heat  exchange 
of  building  surfaces  with  the  sky.  This  apparatus  is 
very  useful  for  the  evaluation  of  nighttime  cooling  of 
low  cost  houses.  It  was  agreed  that  Mr.  Bertolo  would 
send  one  sample  device  to  NBS  for  testing  on  a U.  S. 
building.  Mr.  Bertolo  has  also  obtained  temperature 
transfer  functions  responding  to  a unit  step  change 

of  thermal  load.  What  is  most  needed  for  the  heating/ 
cooling  load  calculation  procedure  of  the  National 
Bureau  of  Standards  computer  program  is  the  inverse  of 
the  temperature  transfer  function;  the  thermal  load 
transfer  function  corresponding  to  a unit  change  of 
room  temperature.  This  transfer  function  is  called  in 
the  United  States  the  weighting  factors.  Appendix  C 
describes  how  the  room  temperature  change  weighting 
factors  were  derived  from  the  data  of  Mr.  Bertolo. 
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6.  Prototype  Desi^^n 


CSTB,  jointly  with  EDF  (Electricite  du  France),  now  has 
design  drawings  for  a prototype  apartment , which  has  been 
conceived  to  embody  the  following  objectives: 

a.  Use  conventional  construction  practice 
bo  Electrically  heated,  but  the  total  con- 
struction cost  of  the  building  is  not  to 
exceed  by  20%  the  cost  of  non-e lectrical ly 
heated  buildings 

c.  In  order  to  justify  the  207o  additional  ex- 
penditure, the  thermal  comfort  throughout 
the  year  is  to  be  better  than  the  conven- 
tional design 

d.  A prototype  building  should  be  completed 
within  two  years 

It  is  believed  that  all  of  these  objectives  can  be  met  by 
a prototype  design  system  which  makes  use  of  a resistance 
heater  imbedded  in  the  floor  slab,  a supplementary  heater 
in  the  ventilation  duct,  exterior  shading  devices,  16  cm 
thick  partition  wall,  and  light  weight  panel  exterior  walls. 
A better  than  usual  indoor  thermal  environment  is  assured 
by  the  design  to  allow  a minimum  temperature  rise  during; 
summer,  minimum  cold  draft  during  winter,  avoidance  of 
cold  spots  in  windows  (use  of  double  glazed  glass),  small 
temperature  stratification  (2  ®C)  and  effective  means  of 
distributing  the  air  to  the  room„  The  prototype  building 
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design  makes  use  of  a heating  control  scheme  such  that 
the  floor  heating  is  done  during  the  night  and  heat 
is  to  be  stored  for  comfort  during  the  day.  Also,  the 
exhaust  heat  recovery  system  is  to  be  effectively  em- 
ployed . 

The  basic  philosophy  of  the  prototype  is  to  design  the 
building  around  the  electric  heating  system  in  such  a 
manner  that  the  maximum  thermal  comfort  and  minimum 
energy  requirements  are  insured  simultaneously.  Mr. 
Borel,  an  exponent  of  the  prototype  design  is  very 
hopeful  that  the  United  States  would  adopt  the  basic 
design  concept  to  the  experimental  housing  of  the  HUD 
Operation  Breakthrough  program. 

7o  Other  Observations  Include: 

a.  Evaporative  cooling  of  a school  building 

b.  Analog  computer  simulation  of  floor  heated 
flats 

c.  Analog  computer  study  of  room  temperature 
response  functions 

8.  Visits  to  Other  Laboratories  Arranged  by  CSTB 

a.  EDF  (Electricite  du  France)  Research 
Laboratories  in  Ecuellis  near  Moret , a 
very  pleasant  one  hour  drive  from  Paris. 

The  author  met  Messrs  Michel  and  Lebault; 
both  of  them  are  on  the  staff  of  the 
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Application  Department  and  are  very  active 
on  the  use  of  computers  for  environmental 
engineering  related  to  electric  heating 
and  air  conditioning.  The  application  de- 
partment, where  environmental  engineering 
studies  are  made,  consists  of  40  research 
personnel  (10  engineers  and  30  technicians) 
and  operates  on  the  annual  budget  of  5 
million  Francs  (one  million  dollars).  They 
have  well  equipped  research  facilities  in 
the  following  areas: 

luminaire 

heat  pumps 

time  constant  for  thermostats 

induction  unit  cooling  with  elec- 
trical heaters 

room  air  temperature  distribution 

heat  pump  heat  recovery  system 
(package) 

office  heat  transfer  affected 
by  adjacent  rooms 

infrared  surface  temperature 
measurement 
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heat  flow  meter  (TNO'' ) application 
for  the  inner  surface  heat  trans- 
fer coefficient  study 

computer  calculation  of  heating/ 
cooling  load  on  the  GE  time 
shared  system  as  well  as  on 
CDC  6600„ 

Energy  calculation  was  done  for  28  days  per  year, 
representing  two  days  per  month,  one  heating  de- 
sign day  and  three  cooling  design  days. 

b.  CEDRIC  (Centre  d* etude  de  Documentation 

et  de  Decherche  pou  I'industrie  du  Chauffage 
du  Conditsonment  d*air  et  branches  annex) 
in  Liege  Belgium  was  visited  on  July  22, 
Accompanying  on  the  trip  were  Messrs 
Bertolo  and  Croiset  of  CSTB.  We  were 
greeted  by  Professor  Bamy,  Messrs  Lebran, 
Nusgens,  Potier  and  Hannoy. 

Mr.  Lebran  has  just  completed  a PhD  thesis 
based  upon  his  work  on  air  movement  within 
a transparent  prototype  room  of  4m  x 4m  x 
2.7m.  He  has  taken  photographs  of  streak 
lines  produced  by  metaldehyde  particles 
following  air  motion  created  by  various 
heating  devices.  His  studies  include  the 
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recording  of  air  velocity  and  temperature  pro- 
file by  thermister  anemometers  and  resistance 
thermometers.  He  used  1 second  exposure  time 
with  the  f/5  camera  opening  and  3X  film.  The 
pictures  were  taken  from  different  positions 
around  the  room  to  traverse  the  flow  field. 

The  partial  photographs  of  the  room  sections 
were  then  pasted  on  18"  x 12"  sheets  in  a 
mosaic  fashion  to  review  the  overall  flow 
pattern . 

Mr.  Lebran  also  has  a 1/3  scale  model 
of  the  prototype  room  filled  with  Freon 
to  simulate  the  convection  flow  pattern. 

In  addition,  the  following  research  activi- 
ties were  observed: 

Electrolytic  analog  studies  on  under- 
ground pipe  heat  transfer  by  Mr, 

Nusgens 

Finite  element  heat  conduction  calcu- 
lation by  Mr.  Fortier 

Thermal  comfort  and  climate  by  Mr. 

Hannay 

c.  On  July  23,  the  research  laboratory  of  COSTIC 
(Comiti  Scientifique  et  Technique  de 
Lindustriedu  Chaffage)  was  visited  by  Mr. 
Galivier,  M,  Anquez  and  the  author.  The 
COSTIC  laboratory  is  operated  by  vari()us 
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organizations  representing  manufacturers 
and  contractors  and  presently  staffed  with 
60  personnel  in  which  twelve  are  engineers. 

We  were  met  by  Messrs  Cadiergues  (Director) 
and  Thin,  and  shown  their  laboratory  facili- 
ties which  included: 

Wind  tunnel  (100  Km/hr  max  speed) 
for  the  wind  around  building 
model  study 

Dehumidification  coil  testing 

Humidity  calibration  chamber 

Testing  facilities  for  air  con- 
ditioners, radiators,  pressure 
relief  valves,  and  heat  ex- 
changers 

We  found  that  they  are  a very  progressive 
user  of  digital  computers » They  have  IBM 
1130  terminals  connected  to  GDC  6600  for 
multitudes  of  HVAC  type  calculations . The 
most  popular  one  of  them  is  for  the  design 
and  selection  of  radiators  with  the  output 
being  complete  line  drawings  for  the  heating 
system  made  by  a plotter.,  The  service  for 
this  kind  of  computer  calculation  is  3.8 
Francs  per  room.  According  to  Mr.  Cadiergues 
there  have  been  260,000  rooms  designed  with 
this  program  during  1969. 
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9.  At  the  conclusion  of  the  author *s  stay  at  CSTB,  Messrs 
Croiset  and  Borel  wanted  to  discuss  the  future  coopera- 
tive program  of  CSTB  and  NBS . The  proposals  were  drafted 
and  approved  by  Mr.  Blachere,  Director  of  CSTB  and  are 
shown  in  Appendix  D. 
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LAYER  L(I)  K(I)  p(I)  C(I)  RES  (I)  DESCRIPTION  OF  LAYERS 

NO.  ft  Btu  per  (hr)(ft)(F)  lb  per  cu  ft  Btu  per  (lb) (F)  (hr)  (sn  ft)  (F)  per  Btu 
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Table  3 


Space  (Room)  No. 

1 

2 

3 

4 5 

6 

7 

8 

9 

10 

11 

Window , Number 

2 

2 

1 

1 

Orientation 

E 

E 

E 

E 

Size  (ft  X ft) 

fS.75 

8.75 

23.5 

23.5 

Solar  Factor 

' x6  ' 

x6  - 

x6  y 

x6 

(CSTB  Definition) 

Ot 

/o 

43 

43 

43 

43 

. . 2 
Floor  area,  ft 

375 

375 

250 

450  925 

875 

875 

725 

400 

625 

625 

Room  volume 

3750 

3750 

2500 

4500  9250 

8750 

8750 

7250 

4000 

6250 

6250 

cu.  ft 


5 

600 


Floor  weight 
factor  Ib/ft 

Space  temperature,  °F 

Occupancy 

Maximum  number  5 

Activity  level 

Btu/hr,  person  600 

Lighting 

Maximum  KW  ^ 

Maximum  W/f t 

Equipment 

Maximum  KW  ^ 

Maximum  W/ft  .75  .75 


80 


75 


0 10  10  5 5 

600  600  600  600 


400  400 


1.5  1.5  4 1.8  3.1  3.5  3.5  2.9  1.6  2.5  2.5 


.7  1.85  1.75  1.75  1.45  0.8  1.25  1.25 


Note:  return  air  picks  up  a part  of  lighting  heat 
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All  the  Walls  (structural  detail  . . Table  ?) 

surface  solar  heat  absorptance  0.35  (structural  detail  .. 
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all  units  in  **foot" 
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Figure  2 Schedules  for  Occupancy,  Lighting  and  Internal  Equipment 
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of  CSTB  Analog  Computation  For  Room  11 
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Figure  4 Sample  Thermal  Load  Calculations  by  the  CSTB  Analog  Computer 


Appendix  A 


July  6 
8:45  a.m. 

9:30  a.m. 

1:00-3:00 

p .m, 

3 :00  p ,m. 
6 :00  p .m. 

July  7.  8 
8:30  a.m. 
9:15  a.m. 

12  :45  p.m 

2 :00  p .m. 

6:00  or  6 

p .m. 


Dr.  Kusuda*s  Stay  at  C.S.T.B. 
(July  5 - July  29,  1970) 


Stay  Organizer : Service  Hygrothermique  et 

Ventilation  (H.T .V . ) 


Left  hotel  by  car  for  C.S.T.B.  Champ s-sur -Marne 
laboratories 

Arrived  at  C.S.T.B. 

Discussion:  Overview  of  H.T.V.  works 

Program  of  the  stay 

Lunch  at  C.S.T.B.  with  Mr.  Blachere  and  H.T.V.  engineers 

Visit  of  the  C.S.T.B.  laboratories 
Returned  to  hotel 


9.  10.  15 


Left  hotel  for  C.S.T.B.  Champ s-sur -Marne  laboratories 
Arrived  at  C.S.T.B. 

Discussions  and  various  works: 

- Papers  sent  by  Dr.  Kusuda 

- Possibilities  of  C.S.T.B.  analog  computer  and  works 
presently  done  with  it 

- Determination  of  convection  coefficients 

- Other  C.S.T.B.  - H.T.V.  and  NBS  - BRD  studies 

Lunch  at  C.S.T.B. 

(On  July  9,  lunch  with  Mr.  Blachere) 

Carrying  on  of  the  discussions  and  works 
30  Returned  to  hotel 
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July  16 


8:30  a.m. 

Left  hotel  for  C.S.T.B.  Champs -sur -Marne  laboratories 

9:15  a.m. 

Arrived  at  C.S.T.B. 

Discussion : 

- Whole  design  of  protot3rpe  building  (how  to  make 
the  thermal  equipment  and  the  type  of  construction 
fit  together) 

12  :45  p .m. 

Lunch  at  C.S.T.B. 

2 :00  p .m. 

Carrying  on  of  the  discussion 

- Prototype  block  of  flats  with  electric  heating, 
studied  by  C.S.T.B.  and  E.D.F.  (French  Electrical 
Authority).  This  study  seems  to  correspond  to  the 
second  phase  of  Breakthrough  operation. 

6:00  or  6 :30 
p • rn  • 

Returned  to  hotel 

July  17 


8:45  a.m. 

Left  hotel  for  Courbevoie  (West  surburb  of  Paris) 

9:15  a.m. 

Arrived  at  Courbevoie 

Visit:  British  Petroleum  Building  (outside  shading 

device ) 

ESSO  Building  (inside  shading  device) 

11:15  a.m. 

Left  Courbevoie  for  Paris 

11 :30  a.m. 

Visit:  Electrical  heating  by  cables  embedded  in  concrete 

floor  (yard)  - Champs -Ely sees 

12  :30  p .m. 

Visit:  Bid  Pereire 

Flat  heated  by  electricity  (double  system:  floor 

heating  and  convectors) 

1:00  p.m. 

Lunch  in  Paris 

2 :30  p .m. 

Left  Paris  for  Valenton 

3 :30  p .m. 

Arrived  in  Valenton 
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July  2 0,  21 
9:00  a.m, 

10:00  a.m, 

1 :00  p .m. 

3:00  or  4:00 

p .m. 

July  22 


July  23 
9:00  a.m. 
10:00  a.m. 

July  24 
8:00  a.m. 
9:15  a.m. 

10:15  a.m. 


Visit:  Hot  water  floor  heating  and  warm  air  with  inlet 

near  ceiling  opposite  to  the  facade 


Left  hotel  for  E.D.F.  Renardieres  laboratories, 
Moret-sur-Loing  (Seine-et -Marne ) 

Arrived  at  E.D.F. 

Discussion  and  visit  of  laboratories  of  Electrical  Appli- 
cations Department  (A.D.E.). 

Contacts:  Mr.  Michel,  Mr.  Lebault 

Lunch  at  Moret-sur-Loing 

Returned  to  Paris 


Trip  to  Liege  (Belgium) 

Visit  of  C.E.D.R.I.C.  laboratory 
Contact:  Mr.  Burnay,  Director 


Left  hotel  for  COSTIC  Saint -Remy-les-Chevreuses  laboratory 
Arrived  at  COSTIC 

Discussion  and  visit  of  laboratory 
Contact:  Mr.  Cadiergues , Director 


Left  hotel  for  C.S.T.B.  Champs -sur-Mame  Laboratories 
Arrived  at  C.S.T.B. 

Discussion  about  conclusion  of  the  stay 

Left  Champ s-sur -Marne  for  C.S.T.B.,  Paris 
Building  - 4,  Avenue  du  Recteur  Poincare 
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1:00  p.m. 

Lunch  with  Mr.  Blachere 

2 :30  p.m. 

Discussion  with  Mr.  Blachere 

4 :30  p .m. 

Visit  of  the  building  (library  ...  ) 

5 :00  p.m. 

Returned  to  hotel 

July  27 
10:00  a.m. 

Met  with  Mr.  Greenberg,  American  Embassy 

2 :00  p.m. 

Returned  to  C.S.T.B. 

Discussion  on  the  future  cooperative  programs 
Conclusion  of  the  exchange  program 

23 


Appendix  B 


Utilization  of  the  C,S.T.B.*s  Analog  Model  to  Improve  the  Digital 

Computer  Program  of  NBS 

prepared  by 

J,  Anquez  and  Bertolo,  CSTB 


Final  Object 

It  seems  to  be  interesting  to  define  from  now  the  final  object  of 
the  work  in  course . 

1.  The  object  of  the  digital  computer  program  worked  out  by 
NBS  is  to  give  the  survey  offices  and  consulting  engineers 
a tool  for  the  calculation  of  the  climatization  plants 
power.  This  program  might  be  improved  owing  to  the  com- 
putations made  by  the  C.S.T.B,  on  his  analog  model  and 
relating  to  the  following  points: 

a.  better  taking  into  account  of  the  building 
inertia  in  the  weighting  factors  calculation 

b.  computations  in  the  case  where  the  inside 
temperature  is  not  steady 

c.  taking  into  account  of  some  non-linearities 
convection  coefficients  on  horizontal  planes 

d.  working  out  mathematical  models  representing 
controlling  systems. 
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2,  If  NBS  had  also  in  view  to  work  out  a computation  program 


for  annual  consumptions,  C.S.T.B.  might  help  NBS. 

a.  on  the  first  hand  in  the  points  detailed  above 

b,  on  the  other  hand  in  utilization  of  climatic  data, 
especially  insolation  and  long-wave  radiation. 


Working  Program 

The  conditioning  loads  calculation  made  with  T.  Kusuda's  data  has 
been  made,  for  the  C.S.T.Bo*s  part,  with  the  whole  of  the  loads  in  the 
air . 

In  T.  Kusuda*s  calculation  the  loads  are  convective  and  radiative. 
Although  the  C.S.T.B.  finds  maximum  load  values  very  near  of  that  cal- 
culated by  T.  Kusuda  with  however  a difference  in  phase;  this  can  be 
explained  by  the  fact  that  the  internal  loads  are  steady  over  several 
hours;  relating  to  the  climatization  loads  this  is  nearly  a steady  state. 

It  would  be  interesting  to  make  again  such  a calculation  for  a 
single  room,  varying  the  inertia  and  the  dimensions  (three  cases  must 
be  adequate  and  for  well  defined  radiative  and  convective  characteristics 
of  the  internal  loads. 

With  the  same  data  it  can  be  considered  that  the  C.S.T.B.  calculate 
weighting  factors  corresponding  to  a variation  of  the  air  temperature 
in  the  room. 
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Then  T.  Kusuda  might  use  these  last  results  in  calculations  that 
might  be  the  following: 

1.  inside  air  temperature  calculation  in  lack  of  climatization 

2 . inside  air  temperature  and  extracted  heat  amount  (from 
which  the  consumption)  for  well  defined  operating  conditions: 

a.  plant  power  capacity  lower  than  the  maximum  load 

b.  plant  operation  stopping  during  some  periods  (week- 
ends ) 

c.  systematic  limitation  of  the  available  power  during 
the  maximum  load  periods,  and  for  given  characteristics 
of  the  equipment  and  of  its  control  (particularly 
proportional  control). 

For  his  part  C.S.ToB.  might  do  the  same  calculations  by  means  of 
the  analog  method  so  that  the  results  may  be  compared. 

The  calculation  of  the  weighting  factors  directly  usable  into  a 
digital  computation  program  can  be  considered  after  the  utilization 
limit  of  a given  weighting  factor  would  have  been  determined,  that  is 
to  say  the  number  of  cases  to  be  calculated  (the  C.S.T.B.  study  into 
the  building  classification  owing  to  their  unit  step  response  might 
be  used  as  a basis). 

Moreover  the  C.S.T.B.  is  in  possession  of  hourly  real  climatic 
data  over  two  years  (grouped  on  punched  cards)  containing  particularly 
the  total  and  diffuse  insolations  on  an  horizontal  plane.  Perhaps 
might  it  be  possible  to  use  them  for  the  calculations. 
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Appendix  C 


Weighting  Factors  for  the  Room  Temperature  Change 

Since  modern  buildings  have  a high  degree  of  dependency  upon  the 
thermal  loads  which  vary  during  the  course  of  a day,  such  as  heat 
generated  by  lighting,  equipment,  and  occupancy  and  heat  gain  by  solar 
radiation,  it  has  been  a trend  to  conduct  the  energy  estimate  calcula- 
tion on  the  hourly  basis. 

In  the  United  States  it  has  been  customary  to  perform  the  hourly 
thermal  load  calculation  on  the  basis  of  constant  room  temperature. 

Actual  thermal  load  imposed  on  the  equipment  (sometimes  it  is  called 
the  heat  extraction),  however,  has  been  observed  to  be  much  less  than 
thermal  load  computed  under  a realistic  condition  when  the  room  tem- 
perature is  allowed  to  change  as  function  of  indoor  and  outdoor  operating 
conditions . 

A simple  yet  dramatic  example  of  this  point  is  illustrated  in 
Figure  1,  which  depicts  the  indoor  temperature  of  the  experimental 
building  of  the  National  Bureau  of  Standards  and  the  controlled  outdoor 
temperature  imposed  upon  the  building.  The  experimental  building  was 
not  heated  or  cooled  during  the  test. 

Although  the  outdoor  temperature  cycle  was  designed  to  produce  a 
large  diurnal  amplitude,  the  indoor  temperature  of  this  building  stayed 
relatively  constant  to  within  +2  °F  of  the  diurnal  average  temperature. 
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If  the  thermal  load  calculation  was  made  on  the  basis  of  constant 


room  temperature,  this  building  would  have  required  heating  during  a 
part  of  the  day  and  cooling  during  another  part  of  the  day.  But  if 
+ 2 °F  change  from  the  average  air  temperature  is  tolerated,  this 
building  would  have  needed  no  energy  for  maintaining  the  required  in- 
door environn^nt. 

A point  to  be  emphasized  is  the  fact  that  if  the  thermal  load  cal- 
culation is  to  be  made  for  the  purpose  of  estimating  hourly  energy  re- 
quirement, it  must  take  into  account  the  effect  of  room  air  temperature 
change . 

Exact  calculation  technique  to  account  for  the  room  air  temperature 
change  upon  the  thermal  load  is,  however,  extremely  complex  and  it  re- 
quires solution  of  simultaneous  equations  describing  the  heat  exchange 

•k 

among  various  surfaces  and  that  between  the  air  and  the  surfaces  . If 
each  complex  calculation  is  to  be  repeated  for  every  hour,  particularly 
for  a building  with  a multitude  of  rooms,  the  computation  time  will  be 
excessive.  In  order  to  decrease  the  computational  effort,  the  ASHRAE 
Task  Group  on  Energy  Requirements  suggests  the  use  of  weighting  factors 
which  may  be  defined  as  follows: 


G,  P,  Mitalas  and  D,  G.  Stephenson,  Room  Thermal  Response  Factors, 
ASHRAE  Transactions,  1967,  Vol,  73, 
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(1) 


where  - effect  on  heating/cooling  load  due  to  the  room  air  tern 

perature  change 


= weighting  factor  for  the  room  air  temperature  change 
= room  air  temperature  change  from  the  set  point 


The  value  of  W,  can  be  determined 
J 


, if  enough  data  for  and  0^  are 


available,  by  a regression  technique,  one  of  which  could  be  the  Wiener 


filtering  technique 

The  data  to  be  used  for  solving  for  from  a regression  procedure 
on  equation  (1),  however,  must  be  obtained  by  exact  solution  of  the 
simultaneous  equation  describing  the  air-surface  heat  exchange,  which 
was  mentioned  earlier.  Such  exact  calculations  have  been  obtained  at 
the  Centre  Scientifique  et  Technique  du  Batiment  (CSTB)  of  France  for 
a typical  room  with  various  construction  characteristics  (Figure  2 shows 
the  room  dimension).  The  exact  solutions  were,  however,  obtained  in 
the  form  of  room  temperature  responsive  to  a step  change  of  thermal  load, 
which  is  schematically  shown  in  Figure  3,  When  the  data  of  CSTB  were 
rearranged  and  plotted  in  semi-logarithmic  sheet  against  time  steps, 
they  form  a linear  relationship  except  for  the  first  value,  such  as 
shown  in  Figures  3,  4,  5,  and  6. 


E.  A,  Robinson,  Multi-channel  Time  Series  Analysis,  Holden-Day,  1967. 
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i 


Using  this  semi- logarithmic  relationship,  the  temperature  response 
function  for  a step  increase  of  the  load  may  be  expressed  as  follows: 


at  t = O,0  = A 

o 

t = A,  0^  = Br 

t = 2A,  0^  = Br 

t = jA,  ^ = Br- 
J 

T<x>  - T. 
where  0.  = ^ 


(2) 


j Tco  - T 


= a set  point  temperature  of  the  room 
Too  = final  temperature  of  the  room  when  a step  of  L = 1 is  applied 

_ -ba 

r - e 

3 = slope  of  the  semi-log  plot  of  the  temperature  response 
A = time  increment 

B = intercept  of  the  semi- logarithmic  plot  at  t = o. 

Table  I shows  the  value  of  A,  B and  r determined  from  the  CSTB  data. 

By  applying  the  equation  (1)  to  the  CSTB  data,  the  following  re- 
lationships are  evident 


1 = W*  0 
o o 


1 = W ' 0,  + W/  0 
o 1 1 o 


1 = W ' 0^  + W/  0,  + W.  * 0 
o 2 1 1 2 o 


(3) 


1 = 2 W . ' 0 , when  W . * = W , (Too  - T ) 

J t-j  J J o 

j=o  J J 
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This  equation  may  be  solved  recursively  to  obtain  the  values  of  as 
follows : 

W*  = 1/A 

0 

W = 7 [1  - W'  . Br]  (4) 

1 A o 

W*  = ~ [1  - W/  • Br  - W*  . Br^] 
z A 1 o 

It  turns  out,  however,  that  the  value  of  decays  very  slowly  as 
the  index  of  j increases  so  that  a large  number  of  weighting  factor  terms 
are  required  to  calculate  by  equation  (1).  The  calculation  can  be 
simplified  very  much  by  applying  a z -transformation  technique  . 

The  z-transformation  presentation  of  equation  (1)  is 


OO  _ C30  , CD 

S Q.  zJ  = ( j;  w.  z^)  ( r.  e.  z^)  (5) 

J •-  J J 

J=o  J-O  J-O 

By  noting  the  fact  that  a polynomial  of  infinite  length  can  be  repre- 
sented by  the  ratio  of  two  finite  length  pol3momials,  one  can  write 

Na 

T.  a.  zJ 

00  . J 

S W. ' Z^  = j=o 

^ iJb  . 

S b.  zJ 
j=o  ^ 


R.  J.  Schwarz  and  3.  Friedland,  Linear  Systems,  McGraw  Hill  Book  COo, 
1965,  pp.  234-267. 
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Simplest  polynomials  to  satisfy  the  CSTB  data,  which  can  be  rep- 
resented by  three  basic  parameters.  A,  B,  and  r,  are  found  to  be  (after 
several  trials) 


s w. ' 
j=0  J 


a 

o 


1 + b. 


+ a^  Z 
Z + b,  z^ 


where  a = 1/A 
o 


r 


b 

o 


-d  + r - j r) 
r(l  - |) 


In  other  words,  the  weighting  factors  for  j = 0,  1,  , , , oo  can  be 
generated  by  these  five  basic  constants  by  a polynomial  division.  The 
use  of  is,  however,  unnecessary  because  the  hourly  load  can  now  be 
obtained  by  a simple  finite  expression  such  that 


Ql  = (a^  - b^  Q^_2)/b^ 


This  expression  requires  temperature  history  of  one  time  step  and  the 
load  history  of  one  and  two  time  steps. 

Tables  II  through  12  show  values  of  a^,  a^,  b^,  b^  and  b^  and 
calculated  by  equation  (8)  for  0^  = 1 and  0^  = 0 for  at  t > o. 
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Also  shown  in  the  same  table  is  the  value  of  0^  calculated  by 


= r + '’i  + ”2  <^t-2  - "1 


(9) 


which  is  the  reciprocal  relation  of  (8)  for  = 1 for  all  t > o. 

Table  II  shows  that  under  a step  increase  of  temperature,  the  value 
of  rapidly  converges  to  a final  stationary  value,  which  is 

1 - r 


Qco  — 


A (1  - r)  + Br 


(10) 


while  the  value  of  9^  slowly  decreases  to  zero  under  the  step  increase 
of  Q. 
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Room  Type 

U 

A 

B 

r 

a 

lo03 

0.9 

0.84 

.841 

b 

II 

0.9 

0.63 

.782 

c 

II 

0.92 

0.90 

.846 

d 

II 

0.92 

0.84 

.739 

a 

2.40 

0.86 

0.85 

.892 

b 

II 

0.86 

0.78 

.790 

c 

II 

0.84 

0.73 

.881 

d 

II 

0.84 

0.50 

.849 

a 

5.88 

0.75 

0o74 

.559 

b 

II 

0.75 

0.71 

.345 

c 

II 

0.71 

0o51 

.646 

d 

II 

0.71 

0.30 

.589 
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Figure  lA  In-  and  Outdoor  Temperature  Profiles  of  a Non-Air- 
Conditioned  Building 


Ceiling 


Room  Air  Change  per  Hour 


N=  1.0 


Figure  2A  Room  Used  by  CSTB  Analog  Computers  for 
the  Calculation  of  Transfer  Functions 
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A 


i Load 


Figure  3A  Temperature  Response  for  a Step  Increase 
of  Thermal  Load 


i 
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Exterior  wall  K = 2.4  watt/m?  ®C 


0.1 1 1 1 1 1— 1 1 1 1 

0 12345678 

f = fime,  hrs 


Figure  4A  Room  Temperature  Response  Function  for  the  inrr-ase  of 

Thenficti  ^ad  Iiiput 
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Exterior  wall  K = 2.4  watt/m?  ®C 

(U=0.424  Btu/hr,  ftf  ®C) 


rigure  5A  Room  Temperature  Response  Function  for  the  Step  Incrcasf 
of  Thermal  Load  Input 
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Exterior  wall  K factor  =1.03  watt/m?  *0 

(U  = 0.I76  Btu/hr,  ft?  ®F) 


t - time,  hrs 


figure  6A  Fvoom  Temperature  Response  Function  for  the  Step  Increase 
of  Thermal  Load  Input 
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Exterior  wall  K factor  = 1,03  watt/m?  ®C 

(U  = 0.176  Btu/hr,  ft?  ®F) 


0 I 23  45  67  8 

/ - time,  hrs 


rc  7A  Room  '1  eiape i ature  Response  Function  for  the  ftep  Incrcasf 
of  Thermal  load  Input 
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Exterior  wall  K = 5.88  watt/m?  ®C 


Figure  8A  Room  Temperature  Response  Function  for  the  Step  Increase 
of  Thermal  Load  Input 
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Exttrior  wall  K = 5.88  watt/mf  *C 

(U  = l.03  Btu/hr,ftf  •F 


Room  Temperature  Response  rin«c  t i>  r for  tSc  Step  Increas 
OL  Thermal  Load  Input 


Figure  9A 
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Appendix  D 


U.  S.  - French  Cooperative  Program  in  Building  Technology 

Proposal  for  Future  Activities  on  Environmental 
Engineering  After  Dr.  Kusuda*s  Stay  at  C.S.T.B. 

During  July,  1970 

1 . Use  of  C.S.T.B.  Analog  Computer  to  Supplement  the  Digital  Computer 

Work  of  N.B.S. 

1.1  C.S.T.B.  would  provide  N.B.S.  the  room  temperature  weighting 
factors  for  the  digital  computer.  The  use  of  the  analog  com- 
puter is  to  expedite  the  calculation  to  account  for  the  room 
heat  capacity  (partitions,  floors,  ...  ) in  the  case  inside 
temperature  is  varying. 

1.2  N.B.S.  would  provide  current  information  on  system  simulation 
in  U.S.A.  (including  characteristics  of  components)  to  C.S.T.B. 
so  that  C.S.T.B.  can  start  a similar  simulation  program  on  their 
analog  computer. 

1.3  C.S.T.B.  would  help  N.B.S.  in  taking  account  of  environments 
control  system  in  the  building  thermal  calculation  by  digital 
computer. 

2 . Climatological  Data 

2.1  C.S.T.B.  would  provide  N.B.S.  one  sample  apparatus  for  measuring 
long  wave  radiation.  As  a start  C.S.T.B.  would  furnish  drawings 
and  schematic  diagrams  of  the  apparatus. 

2.2  N.B.S.  desires  data  for  diffuse  sky  radiation  during  cloudy  days. 
C.S.T.B.  would  explore  possibilities  of  extracting  such  data  for 
their  future  solar  radiation  studies. 
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2.3  N.B.S.  would  provide  to  C.S.T.B.  their  studies  on  earth  tem- 
perature analysis. 

3 . Convection  Inside  Rooms 

Work  is  currently  in  progress  at  N.B.S.  and  is  foreseen  at  C.S.T.B.; 
other  laboratories  (specially  C.E.D.R.I.C.  in  Liege)  have  also  been 
active  in  this  field. 

Everyone  will  keep  himself  informed  and  will  cooperate. 

4 . Ventilation 

N.B.S.  is  interested  in  the  ventilation  processes  used  in  French 
government  sponsored  low  cost  housing  for  the  past  years.  C.S.T.B. 
will  provide  to  N.B.S.  useful  information  in  this  area. 

5 . Overall  Design  of  a Prototype  Apartment  Building  Heated  by  Electricity 
(Building  Designed  Around  Its  Equipment) 

C.S.T.B.  will  provide  N.B.S.  background  material  for  the  general 
ideas  of  the  studies,  the  choices,  specially  in  matter  of  investment 
and  consumption  costs.  N.B.S.  will  see  to  what  extent  the  French 
prototype  design  will  meet  U.  S.  requirements.  N.B.S.  would  explore 
the  feasibility  of  including  the  prototype  concept  for  the  second 
phase  of  Operation  Breakthrough  (O.B.T.).  If  feasible,  a joint 
study  of  N.B.S.  and  C.S.T.B.  could  be  expected  to  adapt  closely  the 
French  prototype  design  to  U.  S.  requirement,  particularly  by  in- 
corporating inexpensive  air  cooling  systems. 
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6 .  Heat  Pumps 


N.B.S.  would  provide  C.S.T.B.  a set  of  technical  literature  per 
taining  to  the  use  of  heat  ptanps . 


N.B.S,  would  help  to  organize  a stay  of  C.S.T.B.  engineers  in  the 
United  States  to  study  heat  pumps  in  dwellings. 

7 . Solar  Factor  Standard 

N.B.S.  would  participate  in  preparing  the  standardized  calculations 
of  Solar  Factors  (S)  planned  by  different  European  countries. 

8.  Trans lat ions 


In  the  frame  of  translation  program,  C.S.T.B.  will  send  to  N.B.S.: 

" two  existing  translations: 

Cahier  n°  334:  "Thermal  weak  points  on  cold  bridges". 

(J,  Berthier)  - translated  by  N.R.C.  Ottawa 

Cahier  n°  608:  "Protection  of  windows  against  solar  radia- 

tion, Explanatory  and  analytical  papers". 

(J,  Anquez,  J.  C.  Borel  & M.  Croiset)  - 
translated  by  B.R.S.  Gars ton 

- three  translations  in  progress: 

Cahier  n°  455:  "The  effect  of  thermal  weak  points  on  the  K 

coefficients  of  sandwich  panels  made  of  con- 
crete and  lightweight  insulation".  (J. 
Berthier  & F.  Clain) 

Cahier  n°  468:  "Protection  of  opaque  panels  by  sun-breakers" 

(J.  C,  Borel) 

Cahier  n°  910:  "Summer  comfort  in  light  structure  school 

buildings",  (J.  C.  Borel),  (livraison  n°  104) 

- two  translations  to  be  done:  Title  V (Thermic)  and  VIII 

(Ventilation)  of  C.S.T.B, *s  Technical  Advice  (Notice  Technique) 
concerning  thermal  problems  and  ventilation. 
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As  soon  as  a new  thermal  study  will  be  published  by  C.S.T.B.  (cold 
bridges,  heat  exchanger,  •..)  it  will  be  sent  to  the  Environmental 
Engineering  Section  of  N.B.S, 

9.  C.S.T.B.  Engineers  Stays  in  the  United  States 

9.1  As  Mr.  Anquez  will  take  part  in  the  "Symposium  on  the  Use  of 
Computers  for  Environmental  Engineering  Related  to  Buildings" 
at  N.B.S.  during  November  30  - December  2,  1970,  it  would  be 
useful  that  Mr.  Bertolo  accompanies  him  and  stays  after  the 
symposium  for  two  or  three  weeks  at  N.B.S.  to  carry  on  the 
collaboration  work  begun  at  C.S.T.B.  in  July,  1970. 

9.2  A short  stay  of  one  or  two  engineers  (two  or  three  weeks) 
devoted  to  heat  pumps  topics  might  occur  in  the  spring  of 
1971;  C.S.T.B.  would  prepare  the  program  with  N.B.S.  help 
in  order  to  firm  up  the  details  at  the  beginning  of  1971. 

10.  N.B.S.  Engineers  Stays  in  France 

N.B.S.  would  send  one  or  two  engineers  for  a period  of  two  or  three 
weeks  to  observe  in  depth  the  work  being  done  in  air  convection  at 
C.E.D.R.I.C.  (Liege),  in  prototype  design  at  C.S.T.B.  (Paris,  in 
physiological  response  at  C.E.B.  (Strasbourg),  in  cold  bridge 
problems  in  prefabricated  walls  at  C.S.T.B.  and  in  thermal  problems 
in  light  structure  school  buildings  at  C.S.T.B. 


Proposed  by: 

Dr.  T.  Kusuda  for  N.B.S. 

Mr.  M.  Croiset  for  C.S.T.B. 

Paris,  July  27,  1970 
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Appendix  E 

Visit  to  Technische  Hogeschool,  Delft  (TNO-TH),  Holland 

The  author  was  met  by  Professor  A.  W.  Boeke  and  Mr.  P.  Euser  of  the 
Mechanical  Engineering  Department  (Werktnigbouwkunde ) ; both  of  them  have 
been  active  in  the  use  of  computers  for  the  environmental  engineering 
calculation.  Prior  to  joining  TH,  Professor  Boeke  had  been  working  with 

Svenske  Flakt  fabriken  (a  Swedish  farm),  where  he  developed  several  of 
the  most  widely  used  European  computer  programs  for  heating  and  cooling 
energy  calculations.  The  details  of  the  program  were  published  in  an 
article  entitled  "New  Development  in  the  Computer-Design  of  Air  Con- 
ditioning Systems"  in  the  October  1967  issue  of  JIHVE  (Journal  for  the 
Institute  of  Heating  and  Ventilating  Engineers).  The  program  is  to  re- 
port a comprehensive  energy  use  for  typical  days  in  a month  for  each  of 
12  months.  The  calculation  for  four  typical  days  were  so  selected  to 
represent  sunny  days  and  cloudy  days  of  weekdays  and  holidays  in  addi- 
tion to  the  design  day  calculation.  The  computation  is  carried  on  the 
module  basis  in  that  each  module  consists  of  many  rooms  with  the  same 
exterior  wall  characteristics.  The  program  permits  the  evaluation  of 
room  temperature  fluctuations  as  well  as  the  room  thermal  response. 
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Mr.  Euser  is  conducting  analog  computer  simulations  of  Phitatron 
(environment  controlled  room  for  the  plant)  and  said  that  the  analog 
system  is  still  useful  for  the  detailed  study  on  control  simulation. 
The  simulation  included  the  simultaneous  transfer  of  heat  and  moisture 
at  the  plant  surface  and  soil  surface.  This  research  is  done  at  TNO 
(Applied  Scientific  Research)  under  the  contract  from  the  Ministry  of 
Agriculture . 
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Appendix  F 

Visit  to  the  Building  Research  Station  (England) 

On  July  31,  the  Environmental  Design  Engineering  Division  of  the 
Building  Research  Station  (England)  was  visited.  The  Division  is  the 
equivalent  of  our  Sensory  Environment  Branch  but  much  larger  in  staff 
(66)  and  activities  (9  sections).  The  author  was  received  by  Messrs 
N.  Milbank  and  P.  Petherbridge  in  the  morning  to  discuss  their  energy 
usage  analysis  and  computer  applications.  They  have  studied  several 
commercial  buildings  with  respect  to  their  weekly  energy  usage  and 
found  that  the  major  portion  of  the  energy  was  spent  for  the  lighting 
and  energy  distribution  (pump  and  fans)  and  that  accounted  for  heating 
and  cooling  was  rather  small.  They  are  currently  continuing  this  survey 
effort  and  have  added  a task  to  include  the  survey  of  expense  for  the 
maintenance  of  equipment  by  sending  standard  questionnaires  to  the  opera- 
tors of  30  buildings.  They  are  finding  that  with  some  scatter  the  main- 
tenance cost  and  the  energy  usage  cost  are  both  4 shillings  per  square 
foot  (707o  of  gross  floor  area).  Their  computer  activities  are,  however, 
just  beginning  and  they  have  recently  completed  a program  to  obtain  solar 
energy  tables.  Mr.  Petherbridge  told  the  author  that  Mr.  Billington  of 
the  Heating,  Ventilating  Research  Association  has  been  active  in  the  use 
of  computers  for  energy  calculation.  Mr.  Billington  has  recently  made 
comparative  analysis  of  Westinghouse,  APEC,  Phillips,  SF  and  Faber  com- 
puter programs  on  a fictitious  building.  The  maximum  thermal  load  cal- 

2 

culated  by  these  programs  varied  widely  (watts/m  of  floor  ranging  from 
3000  to  5000). 
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In  the  afternoon,  the  author  was  met  by  Mr.  Loudon  v7ho  has  been  ac- 
tive in  heating  and  ventilating  research.  Mr.  Loudon  is  currently  studying 
natural  ventilation  of  school  buildings  as  affected  by  outdoor  wind  ve- 
locity. This  study  is  being  conducted  in  a room  (2.53m  x 2.55m  x 2.46m) 
in  a rotatable  experimental  building  with  the  use  of  tracer  gas  (N^O  and 

SF, ).  They  found  that  N„0  is  less  harmful  than  C0„  and  think  that  SF 
o z z o 

should  be  even  better  for  the  study. 

The  natural  convection  study  with  the  use  of  metaldehyde  in  a trans- 
parent room  reported  by  Mr.  Daw  in  their  previous  research  report  was 
terminated  long  ago.  Mr.  Sexton*s  group  is  still  continuing  the  use 
of  the  wind  tunnel  for  the  study  of  air  flow  around  multiplicity  of 
commercial  buildings.  They  are  now  finding  that  the  porosity  of  the 
building  seems  to  alter  the  flow  characteristics.  This  is  interesting 
from  the  standpoint  of  correlating  the  outdoor  air  velocity  profile 
with  the  air  infiltration  characteristics  of  the  building. 
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